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We study the formation of spot patterns seen in a vari- 
ety of bacterial species when the bacteria are subjected to 
oxidative stress due to hazardous byproducts of respiration. 
Our approach consists of coupling the cell density field to 
a chemoattractant concentration as well as to nutrient and 
waste fields. The latter serves as a triggering field for emis- 
sion of chemoattractant. Important elements in the proposed 
model include the propagation of a front of motile bacteria 
radially outward form an initial site, a Turing instability of 
the uniformly dense state and a reduction of motility for cells 
sufficiently far behind the front. The wide variety of patterns 
seen in the experiments is explained as being due the vari- 
ation of the details of the initiation of the chemoattractant 
emission as well as the transition to a non-motile phase. 



Over the past few years, there has been a significant 
increase in our understanding of how spatial patterns 
emerge via the propogating interfacial dynamics of non- 
equilibrium systems [QJ|] . We have discovered how seem- 
ingly different systems can nonetheless exhibit strikingly 
similar behavior, due to the existence of noncquilibrium 
pattern selection principles. These principles rely on the 
idea that the final structure which emerges from an ini- 
tially (linearly) unstable state is affected mostly by the 
nature of the instability, the possible existence of sta- 
ble, highly noncquilibrium isolated steady-state struc- 
tures (such as the single dendrite in nonequilibrium solid- 
ification), and the competition between globally ordererd 
arrangements of these structures as compared to more 
disordered morphologies. This framework has been ap- 
plied to a variety of systems of physical, chemical and 
most recently [p[ biological interest. 

In this workB, we analyze spot and stripe patterns 
seen in bacterial growth experiments; the first such re- 
sults were due to Berg and Budrene || in e. coli. They 
found that cells could aggregate chemotactically, result- 
ing in a wide variety of different colony structures rang- 
ing from arrays of spots to radially oriented stripes to 



arrangements of more complex elongated spots. In their 
study, e. coli were grown on single carbon source me- 
dia and the appearance of patterns only in the case of 
highly oxidized nutrient suggested that respiratory by- 
products (leading to oxidative stress Q) trigger the ob- 
served chemotactic behavior. Two of their figures are 
reproduced in Fig.l. It is worth noting that this class of 
patterns is not limited to e. coli; similar structures have 
been seen in Salmonella typhimurium [Q and in Bacillus 
subtilis We suspect that this class of patterns is a 
universal "possibility" for microbial systems aggregating 
in the face of adversity. 

Normally, bacteria divide and spread out into regions 
of initially low density. This can take the form of an ex- 
panding circle |J, or if some metabolic factor is in short 
supply, a (set of) expanding ring(s) Jl(J. In the pres- 
ence of the aforementioned oxidative stress, the bacteria 
begin to emit a chemoattractant (believed to be aspar- 
tate) which causes them to aggregate via biasing their 
motion. Eventually, the bacteria turn non-motile, freez- 
ing the pattern into place. From the modeling point of 
view, the coupling of the chemoattractant diffusion equa- 
tion to the bacterial density evolution leads to a Turing- 
like instability of the uniform density state. If we were 
to start the entire system at uniform bacterial density, 
patterns would develop as soon as the concentration of 
waste reaches a threshold value. These patterns would be 
defect-ridden, governed by the precise details of the ini- 
tial conditions. This is indeed what was observed in the 
case of purposeful addition of hydrogen peroxide. The 
symmetric structures seen during growth from a single 
inoculated site are thus due to the interplay of the ex- 
pansion with the Turing instability as triggered by the 
waste field, and possibly autocatalytically by the attrac- 
tant field itself. Below, we will see how this occurs within 
our model. 

Based on the above, we propose the following set of 
continuum equations for this system [Tl| ] : 

ft = D p X7 2 p + G(p, n) - v c V • (pVc) - I[w]p 
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p n = I[w]p 

n = DnSI 2 n — apn 
w = D w V 2 w + a w pn — f3 w w — jpw 
c = D c S7 2 c + T[w,c]p- /3 c c (1) 

(2) 

Here p and p n are densities of motile and non-motile bac- 
teria, w is a concentration of respiratory waste products, 
c is the chemoattractant concentration which emitted by 
the bacterium, and n is nutrient which is eaten by bac- 
teria. The functionals T and I are thresholds to be dis- 
cussed below. We have used the freedom to rescale to set 
some coefficients to unity. 

This equation contains a standard growth term which 
we typically take to be of the form G — rp p n/(n + 
no) — o 3 , with p ~ 1 JL2]; this form reflects the nutrient- 
inhibited growth of cells at low food concentration and 
finite reproduction rate which is achieved in a nutrient - 
"rich" limit. At high density of bacteria the growth 
is limited by the nonlinear term — p 3 . Motion is gov- 
erned by diffusion and by a chemotactic term represent- 
ing the response of the bacteria to a gradient in the at- 
tractant. As already mentioned, it will be important 
to incorporate the fact that sufficiently far behind the 
advancing front, the bacteria differentiate into a non- 
motile form. We assume that the transition occurs due 
to the accumulation of effects due to starvation; specifi- 
cally, the transition to the non-motile phase occurs when 
n c = J (no — n)9(no — n) dt exceeds some threshold value 
n tr , i.e. by setting I[w] = SQ(n c — n tr ) jl3| (&{x) is the 
Heaviside function). The rate of waste accumulation is 
proportional to the nutrient consumption term, and we 
assume that the waste decomposes with a rate dependent 
on the bacterial density, as this is in fact the underlying 
reason for the bacteria to aggregate. Finally, the emission 
of chemoattractant c is proportional to the local density 
of bacteria, and it is triggered by local waste field; specif- 
ically, we assume that the chemoattractant is emitted if 
either w > w$ and c > Co or w > wi , where wi > w$ (for- 
mally, T[w, c] — sign(0(w — w<j)Q(c — cq) + Q(w — wi))). 
The difference between w\ and wq represents possible au- 
tocatalytic behavior of the attractant. 

To understand the structure of our model, it is conve- 
nient to first consider the case of constant uniform nu- 
trient (n > no, a = 0) and in the absence of any waste 
effects (i.e. T = 1). The system now has an unstable 
steady-state p = c = 0. If we start with an initial condi- 
tion of localized bacteria density amidst a sea of p = 0, 
the density will spread. In the absence of coupling to c, 
the front would move at a speed at 2*J ' D p r via the usual 
marginal stability criterion [jl4|| . By continuity, the front 
will continue to expand as long as v c is not too large. On 
the other hand, the non-trivial uniform state is given by 
p = %fr 1 c = 1/ (3. The stability of this state to perturba- 
tions with wave vector q is given by the roots of 



(w + D p q 2 + 2r)(uj + D c q 2 + (3 C ) = v c r x/2 q 2 (3) 

It is easy to see that for large enough chemotactic re- 
sponse v c this system has a band of wavevectors with 
purely real and positive growth rates w. This is just 
the Turing instability described above. It does not re- 
quire different diagonal diffusivities D p vs. D c but in- 
stead relies on there being a large cross diffusivity. Com- 
bining these arguments leads one to expect that generi- 
cally there exists an intermediate range of v c for which 
the bacteria density will propagate outward and create 
an expanding Turing-unstable region. Adding the nutri- 
ent field back in does not alter the above conclusions in 
any qualitatively important manner. The importance of 
waste dynamics will be discussed below. 

To simulate the above equations, we used a split-step 
spectral code on a 128x128 lattice; specifically, the lin- 
ear parts of the evolution were treated by FFT and the 
nonlinear pieces by explicit finite differences. Lattice 
anisotropy which comes primarily from the discretizing 
the nonlinear diffusion terms in eq. (Q) was minimized 
by computing them simultaneously on two lattices which 
differ by 45° rotation fUsf , Also, lattice discretization 
effects can lead to spurious (slightly) negative values of 
the density p; this is presumably due to the lattice trying 
to mimic the possibility in the continuum of having pre- 
cisely localized structures due to the nonlinear diffusivity 
Jl6| . To deal with this we used p > 1 (typically ps 1.5) 
and indeed it greatly reduced this spurious effect without 
any qualitative change in the dynamics. We can imme- 
diately verify the above simple features and proceed to 
discuss the formation of complex ordered structures. 

Let us first focus on colony dynamics in the absence of 
a threshold for chemoattractant emission(T(w, c) = 1). 
In Figure 2, we show a snapshot of a typical simulation 
of the model equations. As the interface propagates out- 
wards, it creates a set of concentric rings. The breakup of 
the rings into spots occurs somewhere behind the front, 
dependent on the specific value of v c \ for small enough 
v c , the rings do not break up. Ordered arrangements of 
the spots are hard to achieve, since the interaction be- 
tween different rings is extremely small. Moreover, the 
appearance of a number of rings behind the outer rings 
is not in accord with the Budrene and Berg experiment. 
In most cases the spots appeared right behind the outer 
ring. 

To account for this effect, we impose the threshold for 
chemoattractant emission. The existence of this strong 
nonlincarity makes the rear part of the outer ring ex- 
tremely sensitive to the structure of the waste and chem- 
ical fields. In particular, the emission is triggered prefer- 
entially at certain angular positions as opposed to along 
an entire ring. This then can cure the two previous prob- 
lems - proper thresholding can lead directly to spots and 
previous spots strongly bias the formation of new spots. 
In fig. 3, we show the results of full model simulations, 
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showing significantly improved agreement with experi- 
ment. We now discuss in turn the threshold choices that 
lead to the differing results in figs 3. 

In figure 3a, we see that spots form at positions such 
that they will eventually form radial rows; these rows are 
quite similar to those in the experimental pattern in Fig. 
la. The reason for this is that the concentration of c 
in the ring of currently aggregating spots is nonuniform 
and due to diffusion this inhomogeneity spreads into the 
rear part of the outer ring. This biases the location for 
the new spots so they are created right in front of the 
previous spots. Once the spots form, they are not al- 
lowed to move very far since as nutrient depletes, bacte- 
ria transform into a non-motile phase. If the chemotactic 
response (v c ) is reduced, instead of radially aligned spots 
radial stripes are formed; this too has been seen exper- 
imentally (see Fig. lb and 3b). We note in passing that 
the brightness of spots in experimental images is greater 
for motile bacteria so in all numerical figures we show the 
gray scale plots of (p + 0.5p n ). 

In the other main experimental spot pattern, the spots 
arrange themselves in a manner such that the new spot 
appears roughly between two pre-existing ones at the pre- 
vious ring. By adjusting the values of cq and w\ we can 
produce a similar pattern as well (see fig. 3c). By signif- 
icant increase of the threshold Co we effectively made the 
residual level of chemoattractant irrelevant, and instead 
the waste threshold w\ plays the major role. We should 
note, however, that this structure seems to be missing 
the visually striking spirals which create the sunflower 
impression in the original BB figure (fig la. of their pa- 
per). Whether this difference represents a shortcoming 
of our model which tends to keep spots along circular 
rings, or, as we feel more likely, that a larger system with 
less anisotropy and perhaps a different set of parameters 
could reproduce this visual impression is a question for 
future work. 

To summarize, we have shown how the interplay of 
front propagation and a Turing-type instability can lead 
to spot patterns similar to those observed in bacterial ag- 
gregation. We have concentrated on generic mechanisms 
which should remain important regardless of the details 
of the explicit biological interactions as these details be- 
come clearer. One important conclusion is that the bio- 
logical mechanism for the "turning on" of chemotattrac- 
tant emission is of critical importance in giving rise to the 
observed colonies. Of course, attempts at more quanti- 
tative comparison would require as input the functional 
dependencies of all the various pieces entering into the 
dynamics. We feel that this level of detailed modeling is 
a reasonable undertaking only after a demonstration that 
physics does indeed have a hope of describing what is go- 
ing on without the need to invoke an incredibly complex 
hierarchy of biological mechanisms. This paper has been 
attempt to demonstrate exactly this point and to thereby 
provide motivation for a future quantitative study. 
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FIG. 1. Two examples of patterns of bacteria e. coli in 
experiments by Budrene and Berg [1]: radial alignment of 
spots (a), radially oriented stripes (b) 
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FIG. 2. A snapshot of bacteria density (p + 0.5p n ) at 
t = 18 within the model with D n = 0.4, D c = 0.2, = 0.4, 
a = 0.12, r = 1.8, /3 C = 1, a w = 0.3, /3 m = 0.1, 7 = 0.2, 
no = 0.5, n tr = 0.03, S = 0.5, u c = 6.0,, thresholds for 
chemoattractant disabled, system size is 40, time step 0.02. 
Concentric rings are formed in the wake of the outer ring 

FIG. 3. (a) Same as in Figure 2 but with enabled thresh- 
olds (wo = 0.3, wi — 0.55, Co = 0.04) and v c = 7; radial rows 
of spots are clearly seen (cf. Figure la); (b) Same as in a but 
with weaker chemotaxis (v c = 6), radial stripes similar to Fig- 
ure lb are seen; (c) Pattern of p + 0.5p n within the model (hi) 
with different thresholds (wo = 0.3, Wi = 0.5, Co = 0.5) and 
higher chemotaxis (v c = 9) , other parameters are the same as 
in b. A staggered alignment is seen far from the center. 
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